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Abstract 



I compute the leading color-octet contributions to the process 7 + p — * 
J/ip + 7 (+X) within the non-relativistic QCD (NRQCD) factorization for- 
malism. In the color-singlet model, J/ip + 7 can only be produced when the 
photon interacts through its structure function, while the color-octet mecha- 
nism allows for production of J/ip +7 via direct photon-gluon fusion. Resolved 
photon processes can be easily be distinguished from direct photon processes 
by examining the fraction of the incident photon energy carried away by the 
J/ip in the event. Therefore, this process provides a conclusive test of the 
color-octet mechanism. J/ip + 7 production is particularly sensitive to the 
NRQCD matrix element which figures prominently in the fragmentation pro- 
duction of J/ip at large p± in hadron colliders. I also examine the predictions 
of the color evaporation model (CEM) of quarkonium production and find 
that this process can easily discriminate between the NRQCD factorization 
formalism and the CEM. 
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In recent years, quarkonium production has become the subject of intense theoretical 
and experimental investigation. This has been motivated by the observation of gross dis- 
crepancies between experimental measurements of J/ip and ip' production at CDF [|IJ, and 
theoretical calculations based on the Color-Singlet Model (CSM). Attempts to understand 
this discrepancy focus on new production mechanisms which allow the cc to be produced in 
a color-octet state and evolve nonperturbatively into charmonium. In this paper, I examine 
the photoproduction of associated J/ip + 7 within two different frameworks for compu- 
tation of quarkonium production: the non-relativistic QCD (NRQCD) factorization for- 
malism and the color evaporation model (CEM). In the NRQCD factorization formalism, 
non-perturbative aspects of quarkonium production are organized in an expansion in powers 
of v, the relative velocity of the QQ in quarkonia. For production of S-wave quarkonia, the 
CSM emerges in the non-relativstic limit v — > 0, and color-octet mechanisms are the lead- 
ing order corrections in the v expansion. In the CEM, the hadronization of QQ pairs into 
quarkonia is assumed to be dominated by long distance fluctuations of gluon fields which 
motivates a statistical treatment of color. 

At leading order, the CSM predicts that J/ip + 7 can be produced in photon-proton 
collisions only if the incident photon is resolved, resulting in J j ip that carry a small fraction 
of the incident photon energy. Color-octet mechanisms allow J/ip+j to be produced through 
direct photon-gluon fusion and the resulting J/ip are typically much more energetic. The 
energy spectra of J/ip then directly indicates the relative importance of color-octet and 
color-singlet production mechanisms. Since the CSM, the NRQCD factorization formalism 
and the CEM each differ in the relative importance of the octet and singlet components, the 
J/ip energy spectra predicted by the three models differ in a qualitative way. The leading 
order calculations in this paper will be subject to considerable theoretical uncertainties due 
to higher order QCD corrections. However, for reasons that will be discussed later, I expect 
the shape of the energy spectra to be well predicted by the leading order calculation even 
though the normalization of the total cross section is not. This makes 7+p — > J/ip+j (+ X) 
an excellent process for testing theories of quarkonium production. 

Before I discuss theoretical predictions for 7 + p — > J/ip + 7 (+ X), I will briefly review 
the deficiencies of the CSM and motivate the NRQCD factorization formalism and the CEM. 
A comprehensive review of recent developments in the theory of quarkonia production can 
be found in ref. 0. 

In the CSM ||, quarkonium is viewed as a non-relativistic, color-singlet bound state of 
a QQ with definite angular momentum quantum numbers. Production (or decay) rates for 
a quarkonium state are computed by calculating the production (or decay) of a QQ pair 
which is constrained to be in a color-singlet state and have the same angular momentum 
quantum numbers as the physical quarkonium. The dynamics of the formation of the bound 
state are parametrized by the wavefunction, or its derivatives, evaluated at the origin. 

Such a description is not entirely satisfactory from either a theoretical or an experimental 
point of view. For instance, a calculation of inclusive hadronic decays of P-wave quarkonia 
in the CSM reveals infrared divergences at 0(a^) M. These divergences cannot be sim- 
ply absorbed into the single existing non-perturbative parameter, the wave function at the 
origin. Furthermore, recent calculations of quarkonium production in the CSM fall below 
the observed production of ip' at CDF by a factor of 30 ||. These developments have led 
to a new paradigm for the computation of quarkonia production and decay based on the 
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non-relativistic QCD (NRQCD) factorization formalism of Bodwin, Braaten and Lepage ||. 
In this formalism, the cross section for production of a quarkonium state H is written as: 

a(#) = Eo.<0|G?|0). (1) 

n 

(A similar expression holds for quarkonium decay rates.) The short- distance coefficients 
c n are computable in perturbation theory. The (0|O^|0) are matrix elements of NRQCD 
operators of the form: 

(0|O^i0) = E£<°I4I# (A) + X)(H(X) + X\K n \0). (2) 

X A 

The K n in eqn. (D) is a bilinear in heavy quark fields which creates a QQ pair in a state with 
definite color and angular momentum quantum numbers. These matrix elements describe 
the evolution of the QQ into a final state containing the quarkonium H plus additional 
hadrons (X) which are soft in the quarkonium rest frame. Since this paper will only deal 
with unpolarized production, a sum over polarizations (A) of the quarkonium is included 
in eqn. (Q). Below I will use a shorthand notation in which the matrix elements are given 
as (0^[ g - ) ( 25+1 Lj)). The angular momentum quantum numbers of the QQ produced in the 
short-distance process are given in standard spectroscopic notation, and the subscript refers 
to the color configuration of the QQ: 1 for a color singlet and 8 for a color octet. These 
matrix elements obey simple scaling laws with respect to v, the relative velocity of the QQ. 
This allows for a systematic expansion in the two small parameters relevant for quarkonia: 
a s {2rriQ) and v. For production of S-wave quarkonia, the results of the CSM are recovered 
in the limit v — *■ 0, since the leading matrix element in eqn. (|XJ) can be related to the square 
of the wavefunction at the origin in the CSM. 

This formalism solves the problem of infrared divergences because the infrared diver- 
gences in the calculation of the short- distance coefficients can be systematically factored 
into the matrix elements appearing in eqn. ([]]). Thus a clean separation of long and short- 
distance scales is obtained. Furthermore, new quarkonium production mechanisms are now 
possible, since it is no longer required that the QQ produced in the short- distance process 
have the same color and angular momentum quantum numbers as the quarkonium state. 
Central to the resolution of the ip' anomaly [0] are color-octet processes, in which the QQ is 
produced in a color-octet state in the short-distance process. For S-wave quarkonia produc- 
tion, the color-octet matrix elements are suppressed by v 4 relative to the leading color-singlet 
matrix element. This is not a very large suppression in the case of charmonium since one ex- 
pects v 2 ~ 0.25. Furthermore, in some instances, the short- distance process which produces 
a color-octet QQ can be lower order in the strong coupling than the corresponding color- 
singlet process. In this case, color-octet and color-singlet mechanisms are equally important 
since the color-octet is suppressed by a factor of ~ v A 7r/a s ~ 0(1). 

The factorization theorem for production of quarkonia is a prediction of QCD which has 
been proven at a level of rigor similar to proofs of factorization for the Drell-Yan process in 
hadron collisions. Its utility as a tool for quantitative description of quarkonia production 
depends on the convergence of the expansion in the parameter v. If the expansion in v is 
sufficiently convergent, then only a small number of the non-perturbative matrix elements in 
eqn. ([!]) need to be retained, and the formalism is predictive as well as rigorous. It is possible 
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that for charmonium the velocity expansion is not well behaved, in which case eqn. (|lj) would 
cease to be of practical value. 

A phenomenological alternative to the NRQCD formalism which has recieved renewed 
attention in the wake of the ip' anomaly is the color evaporation model (CEM) ||. The CEM 
is a simple, intuitive model for quarkonia production which relies on a minimal number 
of undetermined parameters. Like the NRQCD factorization formalism, the CEM allows 
for color-octet production, since no constraints are placed on the color configuration of 
the quarks produced in the short-distance process. The CEM differs from the NRQCD 
formalism in its treatment of non-perturbative aspects of quarkonia production. In this 
model, the hadronization of the QQ pair is charcterized by multiple soft gluon emission 
which is not suppressed by any power of v , and so the treatment of color is statistical. The 
sum of cross sections for production of all quarkonium states is related to the differential 
cross section for producing quark- ant i- quark pairs by: 

1 r 2m o da C c. 

&onium — 77 / &Tfl - (d) 

9 J2m c dm 
A similar formula holds for the production of open charm: 

8 f 2m D da C c f , da C c , AS 
cr op en = 77 / dm- h / dm— — . (4) 

9 J2m c dm J2m D dm 

The variable m is the invariant mass of the cc pair. The factor 1/9 represents the probability 
that the cc pair emerges in a color-singlet state. For production of exclusive quarkonium 
states more non-perturbative information is needed. In the CEM, it is assumed that the 
fraction of quarkonium states that emerge as J/ip is independent of the particular production 
process and can be described by a single parameter: 

CJ/-0 = Pj/tpVonium- (5) 

In this paper, I will compute the production cross sections for 7 + p — > + 7 (+ X) in 
the NRQCD factorization formalism and the CEM. Both calculations will be to leading order 
in a s . The absolute normalization of these leading order calculations is not to be trusted to 
a high degree of precision because of large theoretical uncertainties due to ambiguity in the 
choice of the charm quark mass, scale setting ambiguities, and potentially large higher order 
QCD corrections. Therefore measurement of the total cross section will not be a sensitive test 
of quarkonium production theory until the substantial theoretical uncertainties are reduced. 
However, differential distributions are expected to be less sensitive to these uncertainties. 
As an example, the higher order QCD corrections to color-singlet photoproduction of J/if> 
|| increase the total cross section by a factor of 2 but leave the differential cross sections 
da/dz and da/dp± essentially unaffected. Similarly, leading order computations || within 
the CEM show that the model correctly produces the p± dependence of J/ ip production at 
the Tevatron but a K-factor of roughly 2 is needed to achieve correct normalization of the 
total cross section. 

I first discuss computation of J ftp + 7 production cross sections within the NRQCD fac- 
torization formalism. The short- distance coefficients appearing in eqn. ([I]) can be computed 
by matching a perturbative calculation in full QCD with a corresponding perturbative cal- 
culation in NRQCD. The production cross section of a cc pair with relative three momentum 
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q is computed in full QCD and Taylor expanded in powers of q. In this Taylor expansion, 
the four-component Dirac spinors are expressed in terms of non-relativistic two-component 
heavy quark spinors. The NRQCD matrix elements on the right hand side of eqn. ([!]) are 
easily expressed in terms of the two-component heavy quark spinors and powers of q, and 
then the c n appearing in eqn. ([I]) are chosen so that the full QCD calculations and NRQCD 
calculations agree. A more detailed discussion of the matching procedure, along with ex- 
amples of specific cross section calculations, is found in ref. Here I give the parton 
level differential cross sections relevant for resolved and direct photoproduction of J/ip + 7, 
including both color-singlet and color-octet contributions: 



da 



(i + g ^ J/il> + i) 



6Att 2 eto?ajm c 



' s 2 s\ + t 2 t\ + u 2 u\ 



(6) 
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Air 2 



e q aa s 
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't_ + u + 2(2m c ) 2 s \ (O^figi)) 
u t tu I m r 



(8) 



In these formulae, s,t, and u are the Mandelstam variables; s\ = s — 4m^, t\ = t — Am 2 , and 
u\ = u — Am 2 . 

Two variables which are more commonly used in photoproduction are z and p±. The 
variable z is defined to be -Pn /p-y • Pn- The four-momenta of the initial state photon, the 
J/ip, and the proton are denoted p y , p^ and p^, respectively. In a frame in which the proton 
is at rest, z represents the fraction of the photon energy carried away by the J/ip. p±_ is 
the magnitude of the J/ ip three-momentum normal to the beam axis. The differential cross 
section for unresolved photoproduction given above is expressed in terms of these variables 
by making the substitutions: 

_ p\ + (2m c ) 2 (l - z) 
z(l - z) ' 
p 2 ± + (2m c ) 2 (l-z) 



u 



p\ 



1 - z 



(9) 



The total cross section is obtained by convolving the differential cross section with the gluon 
distribution function of the proton, f g / p (x): 



O'tat 



dxdtf g / p (x 



da 

'IE 



dzdp 



2 x fg/v( x ) da 
z(l — z) dt' 



(10) 
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The momentum fraction of the gluon, x, is related to z and p± by: 

„_ s _ pj + (2m c ) 2 (l-z) 

where s 7P is the photon-proton center of mass energy squared. For processes involving a 
resolved photon the formulae in eqn. (|9|) are correct provided z is replaced with z = z/y, 
where y is the fractional momentum carried by the parton coming from the photon. The 
total cross section is then obtained by integrating over the parton distribution functions of 
the proton and photon: 

vtot = J dxdydtfg/ p {x)fg/ y (y)^ = J dydzdp\ X ^ a ^^^ V (12) 

In this s/ySryp. Similar formulae hold for processes in which the initial state 

partons are quarks. 

For calculation within the CEM, the matrix elements for j+g — > c+c+7, g+g — ► c+c+7 



and g + g— >c+c + 7 are found in refs. [ pT] , |T2|| . I fix to be 0.45, which is consistent with 



the results of ref . || . A very similar computation ofe + p— > J/^ + 7 within the CEM can 
be found in ref. I1T31 . 



For both CEM and NRQCD calculations, I use the GRV leading order parton distri- 
butions for the proton JE4|] and the photon ||15| . I use the one loop expression for the 
running coupling with Aq C d = 232 GeV and rif = 3. (This is the definition of a s used 
in the GRV parton distribution functions.) The scale for the strong coupling, /x, is given 
by /i 2 = p\ + (2m c ) 2 , where m c = 1.5 GeV. For numerical values of the NRQCD matrix 
elements I use: 

(O^Si)) = 1.1 GeV 3 
(0^( 3 Si)) = 0.0066 GeV 3 
(0 8 J/ YSo)> = 0.04 GeV 3 
(O^CSo)) = -0.003 GeV 3 (13) 



The matrix element (0(^ (^Si)} is extracted from the leptonic J/ip decay in ref. ||16|| . In ref. 
]T7J, the value of (Og^( 3 Si)} is determined from a fit to charmonium production measured 
at CDF. The values of (O^Y^o)) and (O J s ^( 3 P )} are determined from a study of elec- 
troproduction of J/ip found in ref. p!8[1 . The electroproduction results are consistent with 
linear combinations of (Og^( 1 S'o)) and (Og^( 3 Po)) extracted at CDF [17], low energy fixed 



target experiments [LlJ and photoproduction P0[] . The extractions of color-octet NRQCD 
matrix elements use leading order theoretical calculations, and hence the numerical values 
listed above should be regarded as rough estimates rather than exact results. 

The negative value for (Og^( 3 Po)) may at first sight seem unphysical because the tree 
level definition of this operator is positive definite. This could be taken as a sign that the 
leading order extractions of the NRQCD matrix elements are not be trusted. However, in 



ref. 0] it is argued that loop corrections may actually make the matrix elements of the 
renormalized operator negative. The operator 0^( 3 Pq) recieves a quadratically divergent 
loop correction proportional to the leading color-singlet operator 0(^( 3 Si). Perturbative 
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consistency may suggest neglecting such an effect in a leading order calculation because 
this mixing is formally higher order in a s . However, the mixing is enhanced by a factor 
1/f 4 since the operator 0(^( 3 Si) is lower order in the velocity expansion. In any case, 
photoproduction of J /if) + 7 is not very sensitive to (O^^Pq)), since this matrix element 
only appears in resolved color-octet processes whose contibution is heavily suppressed. 

I compute the total cross section for 7 + p — > J /if) + 7 for the range of energies 

typical of the HERA collider, 30 GeV < yfs < 200 GeV. The cuts z < 0.9 and p± > 1 GeV 
are imposed in order to avoid contamination due to higher twist effects and diffractive 
processes. The theoretical prediction of the NRQCD factorization formalism is shown in 
fig. [IJ. Resolved color-singlet processes give the leading contribution to the total cross section. 
Direct color-octet processes give a contribution that is between 20-50% of the total cross 
section for most values of y/s. Resolved color-octet contributions are negligible. Associated 
J/ip + ^y production is only sensitive to the color-octet matrix element {Og^ ( 3 Si)) and very 
insensitive to (O^ (} S Q )) and (O^ ( 3 P )) . In fig. |, I compare the NRQCD and CEM 
predictions for the total cross section. At most values of y/s the CEM prediction is lower 
than the NRQCD prediction, but never by more than a factor of 2. 

I also compute do /dp\ for the NRQCD factorization formalism and the CEM. The shapes 
of thep_i_ distribution of the resolved color-singlet, direct color-octet and resolved color-octet 
components of the NRQCD predictions are all very similar. The CEM p± distribution differs 
from NRQCD only in the overall magnitude. Since the differences in magnitude are already 
reflected in the total cross section, I do not show the p± distributions. 

More interesting are the distributions in the variable z. In fig. |3], I plot the NRQCD 
predictions for da/dz, at yfs = 100 GeV. Color-singlet processes dominate at z < 0.2 and 
are negligible for z > 0.5, while color-octet processes give a significant enhancement of the 
cross section for z > 0.4. In fig. |], I compare the predictions of the CEM and the NRQCD 
factorization formalism. Since the CEM is dominated by a direct production process, the 
cross section in this model is strongly peaked at high z. Resolved color-singlet production of 
J /if) dominates the NRQCD prediction for the cross section, leading to a substantial fraction 
of J /ip produced at low z, in stark contrast to the CEM where production at low z is very 
small. While the color-octet component of the NRQCD prediction gives an enhancement in 
the large z region, the NRQCD distribution is much flatter here than in the CEM. 

It is clear that the z distribution in associated J /if} + 7 production can be an excellent 
tool for distinguishing between different models of quarkonia production. An observation 
of associated J /if> + 7 events with z > 0.5 would be a strong piece of evidence in favor of 
color-octet mechanisms. Observation of an enhancement of J /if) + 7 events at low z would 
rule out the CEM. Photoproduction of J/if> + 7 also serves as an important new cross check 
of the extraction color-octet matrix elements at the Tevatron, since it is sensitive to the 
matrix element (Og^( 3 Si)), which is most important in fragmentation production of J /if} 
at large p±. Photoproduction of J /if) +7 is complimentary to studies of low energy hadropro- 
duction [Tj|, photoproduction |2(J and electroproduction which provide independent 



measurements of (O^ ? S )) and (O^( 3 P )), but not (0 8 J/ ^( 3 5 1 )). 

After this paper was submitted for publication, two preprints appeared which also dis- 
cussed associated J /if) + 7 production. Cacciari, Kramer, and Greco [^l|] present a detailed 
study of photoproduction as well as electroproduction of J/if> + 7, including color-octet 
contributions. (A preliminary announcement of some of their results is given in a talk by 
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Cacciari and Kramer [p!q| .) Their results are in agreement with the results presented in this 
paper. Early calculations of J/ip + 7 production are reviewed in a recent preprint by Kim 
23| . The calculations of J/ip + 7 production in ep collisions appearing in ref . [^3 are per- 



formed within the color-singlet model. Color-octet contributions may be important for some 
of the distributions presented in this paper. In particular, the photon rapidity distributions 
shown in fig. 9 of ref. p3| were shown to recieve important color-octet corrections in ref. 

I would like to thank Bob Fletcher, Adam Falk and Sean Fleming for many enlightening 
discussions. This work was supported by the National Science Foundation under Grant No. 
PHY-9404057. 
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FIG. 1 - NRQCD predictions for 7+p — > J /ip + r ){+X) total cross section as a function 
of \fs. Dashed, dotted, and dashed-dotted lines indicate the contributions from color- 
singlet mechanisms, resolved color-octet mechanisms and direct color-octet mechanisms, 
respectively. The sum of all production mechanisms is indicated by the solid line. 



FIG. 2 - Comparison of the NRQCD and CEM predictions for the total cross section. 
The solid line shows the NRQCD prediction, the dotted line shows the CEM prediction. 



FIG. 3 - NRQCD predictions for da/dz at y/s = 100 GeV. Dashed, dotted, and 
dashed-dotted lines indicate the contributions from color-singlet mechanisms, resolved 
color-octet mechanisms and direct color-octet mechanisms, respectively. The sum of all 
production mechanisms is indicated by the solid line. 



FIG. 4 - Comparison of NRQCD and CEM predictions for da/dz at ^/s — 100 GeV. 
The solid line shows the NRQCD prediction, the dotted line shows the CEM prediction. 
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FIG. 1. 
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